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We examine two point correlation functions involving the trace of the energy momentum tensor 
in five-dimensional dual gravity theories supported by one or more scalar fields. A prescription for 
determining bulk channel spectral functions is developed. This prescription generalizes previous 
work which centered on one scalar field. As an application of these techniques, we investigate the 
bulk spectral function and corresponding sum rule in the Chamblin-Reall background. We show 
that, when expressed in terms of the beta function, the sum rule for the Chamblin-Reall background 
can be written in a form similar to the sum rule in Yang-Mills theory. 



I. INTRODUCTION 

In recent years, there has been great interest in the 
correspondence between certain strongly coupled gauge 
theories and extra-dimensional gravity theories. This 
duality, originally inspired by the AdS/CFT (Anti de- 
Sitter/Conformal Field Theory) correspondence QUI] has 
led to a set of tools which allow for the computation of 
physical observables in some strongly coupled theories. 
These tools have renewed interest in the possibility of 
a nonperturbative description of QCD and the strongly 
coupled matter created in heavy ion collisions, the QGP 
(quark-gluon plasma). Many excellent reviews on the 
vast literature of this subject are now available, among 
them are [SHS]- 

In order to describe theories which appear in the real 
world, such as QCD, one must move beyond scale in- 
variant theories to nonconformal ones. There are several 
known examples of string theory setups which give rise to 
a dual nonconformal field theory [TDHT5] . An alternative, 
phenomenologically based approach involves construct- 
ing an effective five-dimensional gravity background and 
assuming that it has a dual gauge theory description. 
It is also assumed that all of the AdS/CFT machinery 
still works for this phenomenological setup. The sim- 
plest of such backgrounds involve Einstein gravity, and 
one or more interacting scalar fields. Such models have 
been studied extensively in the literature P3H2I]- To 
date, most phenomenological models have included a sin- 
gle scalar field for simplicity. However, it seems necessary 
to include two scalar fields to properly describe both con- 
formal and chiral symmetry breaking [55J [55] . Fur- 
thermore, many string theoretical setups can be reduced 
to effective five-dimensional backgrounds supported by 
multiple scalar fields [371 HE] 

In this work, we will use such multi-scalar gravity 
theories to study two point functions of the energy- 
momentum tensor in the dual thermal gauge theory. 
Such correlation functions are interesting because they 



allow for the computation of transport coefficients (like 
shear and bulk viscosities) . Transport coefficients can be 
read off from the low energy behavior of the spectral den- 
sity (the imaginary part of these correlation functions). 
These transport coefficients are useful in hydrodynamic 
descriptions of the QGP, and are not otherwise calculable 
perturbatively since the QGP is strongly coupled near the 
phase transition. In addition to AdS/CFT, there have 
been attempts to extract transport coefficients from the 
lattice [55113 1| . Typically one needs to include an ansatz 
for the spectral density in order to perform such an ex- 
tractionPJ In this case, insight gained from gauge/gravity 
duality into the behavior of spectral functions at strong 
coupling could be useful. 

A prescription for determining the bulk spectral den- 
sity in a single-scalar gravity dual was given in |14j . One 
of the purposes of the current work is to generalize this 
prescription to a matter sector which includes multiple 
scalar fields. The results of [T3] employed a coordinate 
system where the scalar field was identified with the ra- 
dial coordinate. While this simplified the analysis, it 
made the generalization to multiple scalar fields diffi- 
cult. In the present work, we will choose our coordinates 
and gauges differently so that generalization to multiple 
scalar fields is straightforward. It should be noted that 
a different method for determining the bulk viscosity in 
multi-scalar gravity dual theories was developed recently 
in [33) . however our method allows for the calculation of 
the spectral function at all values of the frequency w. 

As an example application, we will determine the bulk 
spectral density in a particular single scalar setup called 
the Chamblin-Reall background 34 . This background 
is one of the only known nonconformal setups where ex- 
act results for the bulk viscosity are calculable analyti- 
cally. We emphasize that this model is not particularly 
well suited for QGP phenomenology; it has no conserved 
charge, and also has the peculiar feature of both be- 
ing nonconformal and having a speed of sound which 
is independent of temperature. However, recently it was 
found that the dynamics of a more complicated string 
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1 See 1321 for a critical examination of the methodology of such an 
extraction. 
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theory setup were captured by an effective single scalar 
Chamblin-Reall background 35J. This may be an indi- 
cation that there is a nontrivial connection between the 
Chamblin-Reall background and nonconformal deforma- 
tions of more sophisticated string theory backgrounds. 

While the thermodynamics and transport coefficients 
of this model have already been studied in the litera- 
ture [15 1 120* 1 12 H I5SH55 ] . the full bulk spectral density (at 
nonzero frequency) has not been previously presented, 
and we will use the methods developed in the first part 
of this paper to do so here. (This background does not 
utilize the full power of our methods as it only contains 
one scalar field; however, by working in this model we 
will be able to perform consistency checks with previous 
results.) 

Once we have determined the bulk spectral density, we 
take the opportunity to investigate the bulk sum rule of 
[""""j in the context of this model. We notice an intriguing 
fact that for the Chamblin-Reall background, the shear 
and bulk spectral densities are simply related at all fre- 
quencies. We also demonstrate that the bulk sum rule in 
this model does not agree with the bulk sum rule written 
down in |40j and we explain the reason for the discrep- 
ancy. 

Recently, a thorough study of the shear correlation 
function and sum rule in the Chamblin-Reall model was 
completed in [~*T]. The latter half of the present work 
extends this analysis to the bulk sector. 

The paper is organized as follows. In Sec. [IT] we detail 
the multi-scalar gravity dual in which we are working, 
and set up the relevant perturbations which need to be 
analyzed in order to calculate the two point functions. 
In Sees. Ill and |IV| we generalize the work of [H] by 



with associated spectral densities 



detailing our prescription for calculating the bulk spec- 
tral density, p bulk in a multi-scalar gravity dual. This 
involves the on-shell gravitational action which is given 
in Sec. |III| and the solution of relevant equations for the 



perturbations given in Sec. IV Once we have developed 
our method for calculating the bulk spectral density, we 
apply it to the case of the Chamblin-Reall background 
in Sec. |Vj We provide explicit results for the spectral 
density, Euclidean correlation functions, and we also de- 
rive a bulk sum rule in this model. We summarize our 
main conclusions and discuss some possible directions for 
future work in Sec. |VI| Finally, we provide some of the 
technical details of the calculations as well as some useful 
reference formulae in the appendices. 



II. BACKGROUND AND PERTURBATIONS 

In this work, we will be interested in the retarded cor- 
relation functions of a strongly coupled, nonconformal 
field theory, 

Gg af> (w) = -i f d 4 xe lwt ([T^{x),T a0 {O)])e{t) 1 (I) 
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and similarly for the spectral functions. Here rj^ denotes 
the four dimensional Minkowski metric. Kubo formulas 
relate the low energy limit of the spectral functions to 
the shear and bulk viscosities: rj and £ respectively 
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Often in this work, we will be interested in quantities 
which have their vacuum (zero-temperature) part sub- 
tracted. We use the symbol A to denote this subtraction. 
For example, 



AG R (w) = G R {w) - G R (w) T=Q . 



(7) 



We will compute the correlation functions of interest 
using tools provided by AdS/CFT by working in a dual 
gravitational theory. The dual gravitational theory under 
consideration is a five-dimensional multi-scalar theory, 



S = 




R- M a d^ a 



(8) 



Here, n is related to the five-dimensional Newton's con- 
stant, k = 87rC?5. Summation is implied over the index 
"a" which labels each scalar field. This index runs from 
1 to n with n denoting the total number of scalar fields. 
These scalar fields in the bulk correspond to operators in 
the dual field theory. 

The second term in (J***) is a boundary term (the 
Gibbons-Hawking term). In this term, 7^ is the induced 
metric at the boundary, 9^ v is the second fundamental 
form 



(9) 



the unit vector normal to the boundary is denoted by 
jV m , and V M is the covariant derivative with respect to 
the background metric. 

In order to access two point correlation functions in 
the dual field theory, one must add perturbations to the 
background. 
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The quantities <? M „ and $i are the background fields, and 
and <pi are the perturbations. The metric pertur- 
bations are generally classified under SO (2) rotations in 
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the x±,X2 plane, with the momentum q pointing in the 
X3 direction [HI H21 03] . The bulk mode (which is also 
sometimes called the "sound mode" in the literature) is 
the scalar mode containing perturbations which are in- 
variant under such rotations. 

In this work, we will be interested in correlation func- 
tions at zero spatial momentum q, and thus instead 
of SO(2) rotations, we have the full SO(3) symmetry. 
There are four metric perturbations which fulfill the re- 
quirement of rotational invariance: hoo,h],ho5, and /155. 
Latin indices on the metric components are assumed to 
run over the spatial dimensions i = 1, 2, 3. There are also 
n perturbations, one for each scalar field (p a . 

Not all of these perturbations are physical due to free- 
dom under diffeomorphisms: 



(Pi -Kfi- Vd^ l 



(12) 
(13) 



for any vector tp^. (Unfortunately, this is often called 
"gauge freedom", though it is not the same as gauge 
freedom in field theory). In [T3], only one scalar field 
was considered, and the gauge was chosen so that per- 
turbations ip and h 05 vanished. In this work, we will 
employ the radial gauge (hos = = 0) with perturba- 
tions /ig, h\i and ip a nonvanishing. There are two reasons 
for this choice of gauge over that of [M]. First, doing 
the computation in another gauge allows for consistency 
checks between the two calculations. Secondly, and more 
importantly, the choice of the radial gauge allows one 
to rather easily generalize the computation to multiple 
scalar fields. Such a generalization has not been previ- 
ously presented, and it is one of the central results of the 
present work. 

We define the background and perturbations as 



ds 2 =g tt (z) [1+A(t,z)} dt 2 
+ 9xx{z) [1 + B{t, z)] dx 2 
(pi(t,z) = $i(z) + <pi(t,z). 



■g zz (z)dz 2 



(14) 
(15) 



The vector x spans the usual 3-spatial coordinates. The 
coordinate z labels the "extra" radial dimension; the dual 
four dimensional field theory exists on the boundary at 
z = 0. The background is assumed to be static, and 
all time dependence is contained in the perturbations 
A,B,(pi. 

The radial gauge is initially chosen for the purposes 
of writing down the relevant equations and on-shell ac- 
tion. However, it is desirable to ultimately work in terms 
of gauge invariant perturbations, which are invariant un- 
der infinitesimal diffeomorphisms [501 H31 S3] ■ In [M], the 
more general case of perturbations which depend on both 
space and time was examined. In this more general case, 
it was shown that there are n + 1 gauge invariant vari- 
ables where n is the number of scalar fields in the theory. 
These variables were denoted as Zq and Z^. Here, we 
are interested in the special case where the perturbations 
do not depend on the three spatial coordinates; in this 
case, the variable Zq vanishes. We are left with n gauge 



invariant variables, one for each scalar degree of freedom. 
These variables have the form: 



Z$i(z) = ipi{z) - 



®i{z)g xx {z) 

9xx( z ) 



B(z). 



(16) 



Primes denote derivatives with respect to the radial co- 
ordinate z. It is straightforward to check that these vari- 
ables do not transform under the diffeomorphisms ( 12 ) 
and Jl3l 



At first, it is perhaps surprising that the perturbation 
A (which is dual to the operator T 00 in the gauge the- 
ory) cannot be put into a gauge invariant variable. In 
fact, this behavior is consistent with the expected struc- 
ture of the correlators at zero spatial momentum, as we 
will now show. Because A cannot be put into a gauge 
invariant variable, all correlation functions which involve 
T 00 vanish (up to contact terms) in the limit of zero 
spatial momentum. The vanishing of spectral functions 
involving T 00 is expected at q = on general grounds. A 
general correlation function of components of the energy 
momentum tensor can be defined to satisfy the Ward 
identity [42]l43] 



q ll G liVaP = contact terms. 



(17) 



The contact terms which appear on the right hand side 
of this equation depend on how the correlation function 
G is defined. (Note that G here does not necessarily co- 
incide with the retarded correlation function Gr defined 
in ([lj). Here, the only necessary detail concerning the 
contact terms is that they are all real, so that regardless 
of the definition of the correlation function, the spectral 
density is the same as that given in (2]). The Ward iden- 
tity implies that the spectral density is transverse: 



g„ [ Im G^] = -q^pT^ = 0. 



(18) 



Writing this in components and using the definition of 
the four-momentum q^ — (w,q), 



(19) 



Taking the limit where q — > 0, the right side vanishes be- 
cause the correlation functions are nonsingular functions 
at vanishing spatial momentum^] Then, we see that 



p Ouaf) (w,q^0)=0. 



(20) 



Our calculations in the dual gravity theory will be con- 
sistent with this result. 



III. ON-SHELL ACTION 

The prescription for calculating Minkowski space two 
point functions of a strongly coupled field theory from 



2 This is true except in some very special cases where excitations 
at w 5^ do not decay as a function of time 1451 146| 
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gauge/gravity duality was first given in [17]. One needs 
to solve the linearized Einstein equations for the pertur- 
bations, and plug the results back into the action which 
has been expended to quadratic order in the perturba- 
tions. In this section we will expand the action and write 
it in terms of the gauge invariant variables. The main 
result of this section is an expression for the spectral 
densities in terms of the gauge invariant variables. 

Let us assume that we can write our background metric 
in the coordinate system 



ds 2 = g xx {z) 



-f(z)dt 2 + dx 2 + 



dz 2 
W) 



(21) 



This coordinate system is only chosen as a calculational 
aid for the moment. We shall see that our final results 
will be valid in any coordinate system. With the use of 
the background equations of motion, and the linearized 
equations for the perturbations one can write the on-shell 
action, expanded to quadratic order in the perturbations, 
in the form 



S 2 



V3 f dw 



Q 3/2 f 



4k ./ 2tt 

+ Q{w,z) T Z(z)Q(-w 7 z) 



Zq, a (-w,z)Z'$, a (w 1 z) 



(22) 



A derivation of this equation is presented in Appendix 
|B] The perturbations depend on w rather than t now 
because we are working with the Fourier transform. A 
summation over the repeated index a (which takes values 
from 1 to n) is implied. Here, V3 denotes the 3-volume 
(integration over the three spatial coordinates), Q is a 
vector made up of perturbations, 



Q(w,z) 



( A(w,z) \ 
B(w,z) 
<pi(w, z) 

\ (p n (w,z) ) 



(23) 



and £ is a symmetric (2+n) x (2+n) matrix, the compo- 
nents of which are unimportant for our present purposes 
(but are given in Appendix [B]) . This expression for the 
on-shell action does not contain possible counter terms 
that are necessary to regulate UV divergences. We do 
not need to consider such counter terms, because the UV 
divergences generally only affect the real part of the on- 
shell action [TU 27J HSJ. In this work, we will only be 
interested in the spectral functions, which are computed 
from the imaginary part of the on-shell action. In other 
words, regularizing counter terms will change the compo- 
nents of the matrix £, but this matrix does not contribute 
to the spectral functions. 



The imaginary part of the on-shell action can be writ- 
ten 



Im <Sq = 



TsV 3 
4i 



dw 
2^ 



9xx{z) 



■A/2 



/(*) 



9xx{Zh)J f'(zh) 



x Z' 9a (w, z)Z^ a (-w, z)-(w-> -w) . (24) 
Here the temperature T and entropy density s have been 
used (details can be found in Appendix [B]). It is well 
known that the imaginary part of the on-shell action is 
proportional to a conserved quantity (see [TH, |4"T1 EE]), 
and thus we can evaluate the above expression at any 
value of z. We choose to evaluate it at the horizon, and 
make use of the incoming wave condition 



Z' m (z h ){z - z h ) 



-I'll! 



Z^{z h )[l + 0{z-z h )\, (25) 



4ttT 

and the near-horizon behavior of the metric 
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The result is 



Im S 2 = ~ I W -^- [Z<s, a {w, Zh)Zq, a {—w, Zh)] , 



which can be written 



— Im £2 

n T/ 
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(28) 



v dw 
~2tT 



Z% a {-w) 



4W4H 
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We have chosen to write it in this form because the quan- 
tity in square brackets is most directly related to the 
spectral density. The superscripts "0" and "h" denote 
the boundary and horizon value of the perturbations at 
z — > and z — > Zh- Often the boundary values are diver- 
gent (though the final answer will not be). To be precise, 
one should introduce a UV regulator e 



Z% a = Z$ a (w,e), 
Z% a = Z$ a (w, Zh), 



(29) 
(30) 



which is assumed to be small. At the end of the calcula- 
tion one should take e — s- 0. 

We will need to use the fact that due to the definition 



of Zm (16) 
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ZU-w)ZUw) = \Q\-w)] 
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h(e) 2 

-fci(e) 

V 



2 + i 




2 + n 

\ 




^(w) (No sum over i), (31) 



/ 



where 



the spectral density matrix is 



ki(z) 



g'xx(z) 



(32) 



We have added column labels above the matrix to aid the 
reader. The matrix is square, so the same labels apply 
to the rows. To be clear, all ellipses (...) in the above 
matrix denote zeros; the nonvanishing elements of this 
matrix are written explicitly here. 

The Son and Starinets prescription for the correlation 
functions [27] , (see also [15] for a nice treatment of mixed 
operators) can be stated as follows. After writing the on- 
shell action in the form 



Im <S 2 = 



dw 
2tt 



[Q°(-w)} T MQ°(w), (33) 



p{w) = -2M. 



For simplicity, we define 



Ri{w,e) 



Z^i-w^iw) 



(34) 



(35) 



The repeated index i is not summed, and the second 
equality is due to the fact that Z{w)* = Z(—w). Com- 
bining (28 1, (31 1, (34), and (35) we have the main result 



of this section: 
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The ellipses in this matrix do not necessarily denote ze- 
ros. To clarify, the entries in the second row and the 
second column of this matrix are all nonvanishing except 
for the zero which is written explicitly above. Similarly, 
the entries on the diagonal are all nonvanishing, except 
the upper left corner. All other entries are zero. The 
column labels are to aid the reader in distinguishing the 
significance of each entry; the matrix is square so the 
same labels apply to the rows. 

This matrix contains all information about spectral 
densities of the relevant operators in the dual field the- 
ory. For example, suppose one is interested in the spec- 
tral density of the two point function (0 B 0' il ) for two 
operators O b and O* 1 which are dual to the fields B 
and $i. The examination of the above matrix tells us 
the answer is 

p B ^(w) = - lim S ^k x {e)R l {w,e). (37) 



Note that all spectral functions involving the operator 
Oa ~ T 0Q vanish. This is the expected behavior at zero 
spatial momentum as explained in Sec. [n] 

In this work, our primary interest is p bulk . The op- 
erator dual to the perturbation B is \T\ (the factor of 
a half introduces an extra factor of 4 in the correlation 
function) |42| . Thus, we need the matrix element in the 
second row, second column: 

n 

p bu!k (w) = SW y, ka{e) 2 Ra{w ^ £) (38) 

7T c— >0 * — ' 

a— 1 

If one is only interested in the bulk viscosity, one can 
apply the Kubo formula ([6]), 

n 
a— 1 

It is worth emphasizing that the k t functions do not 
change under redefinition of the radial coordinate. Cer- 



G 



tainly the bulk spectral density and bulk viscosity which 
are physical quantities must not depend on the choice of 
coordinates in the gravity dual. Thus, the R a functions 
must also be coordinate independent (we will be able to 
see this explicitly in the next section). Therefore, this 
prescription for calculating the spectral density and bulk 
viscosity is valid in any background coordinate system. 

In order to evaluate these quantities one must compute 
the Ri functions by solving the linearized equations for 
the gauge invariant fluctuations Z^i . In the next section 
we will write down the relevant equations. 



IV. LINEARIZED EQUATIONS 

In the previous section, we gave results for the bulk 
spectral function (and bulk viscosity) in terms of gauge 
invariant perturbations. In order to get an explicit result, 
one must solve the linearized Einstein equations for the 
perturbations. We will now detail the equations which 
need to be solved. 



perturbations is set to zero, the gauge invariant variable 
Zq vanishes identically. Hence, this equation simplifies 
to 



1 



d z [y/^ggr'Z'n] ~ w 2 g tt Z m 




(43) 



d z [^99 zz Kki]}=0. 



There are n such equations, one for each of the Z$i 
variables. This system of n coupled equations must be 
solved in order to determine the bulk spectral density in 
a generic multi-scalar gravity dual. (Usually, the com- 
plexity of the equations requires a numerical solution). 
If one is only interested in the bulk viscosity, one may 
set w — in the above equation. 

The equation ( 43 ) comes from [20] . It is clear from the 



derivation in that paper, that no special coordinate sys- 
tem is assumed other than the usual black brane ansatz: 



ds 2 = gtt(z)dt 2 + g xx (z)dx 2 + g zx {z)dz i 



(44) 



A. Theories with multiple scalars 

In (20j , sound mode perturbations were examined for a 
general, multi-scalar gravity dual with an arbitrary num- 
ber of spatial dimensions denoted by p. The equations 
for the gauge invariant variables were given in full gener- 
ality, assuming both space and time dependence. In the 
general case there is another gauge invariant variable Zq 
which appears in addition to the Z$i which have already 
been introduced in the previous section. The equation 
that is relevant for our purposes i^] 



This justifies our discussion in the previous section; equa- 
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Z d z [V^gg zz Z^] - Z^g zz (w 2 g u 

g 




q 2 9 xx 



w 

7 



(40) 



Here, the notation T>l denotes the logarithmic derivative 



V L [Y(z)]=Y'(z)/Y(z), 
and the quantity a is defined as 



a(z) — q 2 



(p-l) + 



Vl[9u\ 

'Dl [g X x. 



pw 



(41) 



(42) 



In the current work, we will limit ourselves to the case 
of zero spatial momentum q = 0, and three spatial di- 
mensions p — 3. When the spatial momentum of the 



3 This is Eq. (66) in 1201 . Readers should note a different definition 
of the function / in this paper. We have re-written the equation 
to correspond with the definitions given in the present work. 



tion ( 43 1 is valid in any coordinate system with these 



symmetries, and so of course the set of solutions (Ri func- 
tions) is also independent of the coordinate system. 

Let us now summarize the procedure to determine the 
bulk spectral density in this general multi-scalar model. 

1. For the potential V(&i...& n ) under consideration, 
one must solve the background equations of motion 
to determine g xx , g t t, g zz , $i...$ n . Once these are 
known the fc, functions are known. The relevant 
background equations are written in Appendix [Aj 

2. Once the background is determined, one must solve 
the set of equations (|43l subject to the standard 



incoming wave condition. This is usually done by 
inserting the incoming wave ansatz 



Z m (z) = f(z)- iw ^ T Yi(z) 



(45) 



with the assumption that Yi is regular at the hori- 
zon. A convenient way to do this numerically is to 
begin numerical integration at the horizon and inte- 
grate toward the boundary to determine the values 
of Z\^ We refer the reader to jH] as the numeri- 
cal procedure outlined there could be implemented 
here as well. 

3. Once the boundary values Z^ are determined nu- 
merically, the Ri functions are known, and hence 



one can use the formula ( 38 ) to determine the spec- 
tral density. 

This completes our generalization of the work of [2] to 
multiple scalar fields. 



7 



B. Single scalar theories 



A. Bulk spectral density and Euclidean correlators 



If only one scalar field is present, it is possible to com- 
pletely remove the potential from the gauge invariant 
equations, since we can trade derivatives with respect to 
<f> for derivatives with respect to z using the chain rule. 
After doing so and making judicious use of the back- 
ground equations of motion (for details, see Appendix 
[C| , we find the single gauge invariant equation reduces 
to the following form. In the case of one scalar field, 
there is only one gauge invariant variable, and so we have 
dropped all the i subscripts for simplicity: 



1 



The associated spectral density is 

p hulk (w) = —\imk(e) 2 R(w,e). 

7T e-+0 



. (46) 
(47) 



CHAMBLIN-REALL BACKGROUND 



Despite the relative simplicity of ( 46 ) , it docs not seem 



to be possible to solve the equation analytically in gen- 
eral, even in the limit of vanishing w. However, we are 
aware of one special case where certain analytical results 
are possible. This case is referred to as the Chamblin- 
Reall background, where the scalar potential is a pure 
exponential. The potential is defined as 



vm = - 



6 (2-6) 
L2 (1-26^ 



exp 




(48) 



with the numerical factors chosen for later convenience. 
The parameter L has dimensions of length; in the con- 
formal case of 5 = it is the AdS curvature radius. The 
metric and background field are written as follows. 

ds 2 = b 2 (z) \-f(z)dt 2 + d X 2 ' /: ~ 



b(z) = 



1/(1-25) 



(49) 
(50) 



m 
*(*) 



= 1 - 



2(2-<5)/(l-2<5) 



(51) 
(52) 



z 

^125 \og[b(z) 

The conformal symmetry breaking parameter S ranges 
from < 6 < 1/2. It is related to the trace anomaly and 
speed of sound as [301 13Z] , 

S= ^^^M (53) 
2e 2 K ' 

where e, P, and v s denote the energy density, pressure, 
and speed of sound of the dual fluid. Note that the 
Chamblin-Reall background is special, because the func- 
tion k(z) is independent of z. 



In this background, the equation of motion for Z$ be- 
comes 



1 



:d z y^gg zz Z^\ - Z^g zz g tt w 2 = 0. (54) 



V-99 2 

At this point, we notice an interesting coincidence be- 
tween this equation for the bulk perturbation, and the 
case of the shear perturbation. In [5T], the shear spec- 
tral density was examined in this model. In complete 
analogy with the expressions above, it was found that 



shear 



sw 
An 



H(w,z h )H(-w,z h ) 



H(-w,e)H(w,e) 
with H being the solution to the equation 
1 



-d z 



-99 



99 



w 2 g zz g tt H 



0. 



(55) 



(56) 



It is clear that the equations for H and for Z$, are identi- 
cal in this background only. (This fact was also recently 
pointed out in [49 J. Using the fact that H = Z®, we 
find that the bulk and shear spectral densities are very 
simply related; their ratio is a constant: 



„bulk 



, = 4 lim k(zY = 125. 



(57) 



Applying the Kubo formulas to the case at hand, we see 
that 



1 



(58) 



which is a well known result (for example, see [T4l 120) ). 
One of the novel observations of this work is that the 
above ratio between £ and r\ is actually a special case of 
the more general fact that the full shear and bulk spec- 
tral densities are simply related at all values of it;. In 
Appendix |Dj we re-derive this result using the methods 
of [13] and show that both methods give the same answer. 

The shear spectral densities xwere computed in |41) . 
and hence we can use the results given there to compute 
the bulk spectral densities. Some results are plotted in 
Fig. [1} 

The spectral density is related to the full, Euclidean 
correlation function through an integral transform, 



G E (r) = 



1 



dwp(w) 



cosh [w(t - 13/2)} 
sinh[w;,5/2] 



(59) 



where r is the Euclidean time variable, which has period 
j3 = 1/T. Euclidean correlation functions can be com- 
puted on the lattice. In Fig. [2] we plot this quantity 
(with the zero temperature part subtracted) in the bulk 
channel for the Chamblin-Reall background. Despite the 
oscillations in the spectral function, the Euclidean corre- 
lation functions turn out to be smooth. 



Zero-temperature subtracted bulk spectral density 






= 0.10 




= 0.15 


s 


= 0.20 



12 3 4 

w/(2 n T) 

FIG. 1. Plots of the zero-temperature subtracted spectral 
density versus frequency for several values of 5. For large w, 
the spectral density always approaches the zero-temperature 
result. Note that unlike the shear spectral density, the inter- 
cept at w — increases with 5. This is a demonstration of the 
fact that the bulk viscosity increases with 5. In a conformal 
theory, (5 = 0), the spectral density vanishes identically. 



Zero-temperature subtracted Euclidean correlators 
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5 = 0.20 
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FIG. 2. Plot of the zero-temperature subtracted Euclidean 
correlation function associated with the bulk spectral density 
as a function of the Euclidean time r for various values of 5. 
In this model the value of 5 (and hence the value of the bulk 
viscosity in this model) has a strong effect on the shape and 
magnitude of these functions. In a conformal theory, (5 = 0), 
this function is identically zero. 



Multiplying both sides by 12(5, and using (57), we find 
the result for the bulk channel 



36(5 
10-85 



(e + P)= r - [p™\w) - p*t k M] ■ (61) 
Jo w 



In [JO] , the left side of the above sum rule was derived for 
Yang-Mills theories, and was found to take on the value 



3(e + P)(l-3^)-4(e-3P). 



(62) 



If we naively evaluate this quantity in the Chamblin-Reall 
metric using ( 53 1 , we find 



3(e + P)(l-3i>2)-4(£-3P) 



-6<P 
2-5 



(e + P). (63) 



Clearly, there are differences between the bulk sum rule of 
[40] and the bulk sum rule for the Chamblin-Reall metric 
(both the sign and the leading power of 5 are different) . 

What is the nature of this difference? The derivation 
of the sum rule [4QJ employs asymptotic freedom, and so 
we should examine the beta function in the Chamblin- 
Reall background to see whether this model meets the 
necessary requirements. 



C. Beta function and trace anomaly 

One can show that in single scalar gravity duals, the 
trace anomaly is consistent with that of QCD, 



(Tr(P 2 )) =4(e-3P) 



provided that we identify [18] 



(Tr(P 2 )) = X t ^ (O.) 

ff(A t ) = /3 d* 
A t V 8 db' 



(64) 

(65) 
(66) 



Here A t = g 2 N c is the t'Hooft coupling, F 2 is the square 
of the field strength tensor F^ v , and b is defined from 
the metric as written in the coordinate system ( 49 ) . The 
field theory operator dual to the scalar field * is denoted 
as 0$. The t'Hooft coupling is usually identified with a 
scaled version of the dilaton 



Bulk sum rule 



With the spectral density in hand, it is a simple matter 
to investigate the bulk sum rule of [40]. The shear sum 
rule in this background was already studied in [41] . with 
the result 



10-8(5 



(e + P) 



dw 



[p sh ™(w)- P f^(w)} . (60) 



co A t 



(67) 



with Co being an unknown proportionality constant which 
will not enter our equations. Note that 



/?(At) 
At 



d* dz 
dz db 



3*'(- >//. 



>(*) 



9xxi Z ) 



k(z). 
(68) 

In order to get an explicit result for (3 in the general case, 
one would have to invert Af(z) and substitute it above to 
get k(z(Xt))- In the case of the Chamblin-Reall metric, 
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this is unnecessary because k is a pure constant. In this 
background, then 



/?(A t ) = -3A t i 



(69) 



The beta function is negative, and hence this model pos- 
sesses asymptotic freedom. There must be another rea- 
son why the sum rule in this model does not agree with 
that of gD]. 

The left side of the sum rule is, more generally, 



AG^ ulk ( 
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AG^ ulk (w = 0), 



(70) 



where we use the symbol A to denote a zero-temperature 
subtraction as in 0. The low energy piece (at w = 0) 
provides a universal contribution from hydrodynamics, 
so this is not the source of the difference. The crucial 
difference is the source of conformal symmetry breaking 
in each model. In gauge theories, the conformal sym- 
metry is broken by the running of the coupling. If the 
theory is asymptotically free, at large energy scales the 
coupling vanishes due to asymptotic freedom, and hence 
the high energy term (at w — ioo) vanishes. In other 
words, conformal symmetry is restored at high energies 
or temperatures in such gauge theories. 

In the Chamblin-Reall theory, there is a hard breaking 
of conformal symmetry due to a parameter in the La- 
grangian. This parameter does not run with energy scale 
and hence at large energies or temperatures, conformal 
symmetry is not restored, and thus AGr(w — ioo) does 
not vanish in general. Thus, the different form of the left 
side of the sum rule in this theory is a consequence of the 
fact that the conformal symmetry breaking is present for 
all temperatures. 

Despite the fact that the left side of the sum rule ap- 
pears different in each theory, we can find some similarity 
by writing the left side in terms of the beta function. In 
what follows, we will assume that the conformal symme- 
try breaking is small (e — 3P) <C e, P. 



1. Yang-Mills theory 

For gauge theories, the region of small conformal sym- 
metry breaking occurs at weak coupling. The pressure 
can be computed as 



P(T) = T i( A + Bg i + Oig 3 )) 



(71) 



using s = P'(T), e + P = Ts and v 2 = P'(T)/e'{T), one 
finds that, to lowest order in the coupling constant g, the 
left hand side of the Romastschke-Son sum rule is 



3(e + P)(l - 3v 2 ) - 4(e - 3P) » BT l 



dP(g 2 ) 
dT : 



Where 



/V) = T 



dT 



(72) 



(73) 



The beta function is computed as 

P(g 2 ) = -bog* + b ig 6 + 

Again, to leading order, 



T 



dT 



2b Q g 2 p(g 2 ) = 2 



V 9 



(74) 



(75) 



In all, then the left side of the sum rule becomes (to lead- 
ing order in the conformal symmetry breaking parameter, 
<?), 

3(e + P)(l - 3v 2 ) 4( £ - 3P) « ^P (76) 

where Pq denotes the pressure in the conformal limit 
Po = AT 4 . Inserting factors of N c and using the fact 
that the t'Hooft coupling is A t = g 2 N c , we find 



3(e+P)(l-3u^)-4(e-3P) 



2PA t f/3(X t ) 

an; Po {— 



(77) 



2. Chamblin-Reall background 

For the Chamblin-Reall background, the conformal 
symmetry breaking parameter is S. The left side of the 
sum rule is to lowest order: 



18 



Sx4P + O(S 2 ), 



(78) 



Here, again, P denotes the pressure in the conformal 
limit P = P(S — > 0), and we have used the fact that 
e = 3P in the conformal limit. Employing the use of the 
beta function 



/3(A t ) 



2 ' 



(79) 



we find the left side of the sum rule can be written 



16 p fP(\t) 



(80) 



Comparing ( 80 1 to ( 77 ) , we note that when the conformal 



symmetry breaking is small, the bulk sum rule is well 
approximated in both theories as 



CPr 



V A, 



dw 



trH-fiw]. (si) 



This equation is approximate because we have only in- 
cluded the leading order term in the deviation from con- 
formality on the left side. While the functional form is 
similar in both theories, C is theory dependent. For the 
Chamblin-Reall background we have shown C = 16/5; 



for pure Yang-Mills theory [50], it is C 



-5\ t /8ir 2 . 



It would be interesting to see whether this dependence 
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holds in other theories beyond those considered here. It 
is interesting that C contains a different dependence on 
the coupling in each theory; this is a reflection of the 
fact that to lowest order the beta function goes as At in 
Chamblin-Reall, but as A 2 in Yang-Mills theory. 

In the theories we have considered, the left hand side 
of the bulk sum rule has the same dependence (on the 
conformal symmetry breaking parameter) as the bulk vis- 
cosity. In fact, it has been argued on general grounds that 
the bulk viscosity should be proportional to two factors 
of the beta function in [51] . Despite the fact that the 
sum rule (and bulk viscosity) both depend on two pow- 
ers of the beta function in both theories, the dependence 
on (e — 2>P) is different. By looking at the Kubo formula 
for the bulk viscosity one might expect the dependence 
C ~ (e— 3P) 2 . Here, we see this is not generally true. The 
quantity e — 3P is not necessarily linearly proportional 
to the beta function. 



VI. CONCLUSION 



both Yang-Mills theory and in the Chamblin-Reall back- 
ground is proportional to two factors of the beta function. 

When one attempts to extract transport coefficients 
from the lattice, one typically needs to make an ansatz 
for the spectral function. The results from AdS/CFT 
computations like the one presented here could be used 
in this process (though, as we mentioned previously, it is 
necessary to perform the computation in a theory which 
shares more phenomenological features with QCD). 

There are many possible directions for future work. It 
would be useful to further generalize our prescription by 
including other matter fields beyond the simple scalars 
included here. The inclusion of vector fields allows one 
to study holographic models at nonzero chemical poten- 
tial. One could also try to generalize our methods to an 
arbitrary number of dimensions. Furthermore, it would 
be interesting to compute the correlation functions at 
nonvanishing spatial momentum. In this way, one could 
have a complete generalization of the known conformal 
results of [52] ■ We hope to address some of these issues 
in future publications. 



The main results of this work were presented in two 
parts. First, we presented a prescription for calculating 
the bulk spectral density (and hence the bulk viscosity) 
for a general, five-dimensional multi-scalar gravity dual 
theory. In this way, we have generalized the results of [14] 
to multiple scalar fields. In so doing, we have calculated 
the on-shell gravitational action relevant for the compu- 
tation of bulk mode two point correlation functions. We 
have also given explicitly the set of gauge invariant equa- 
tions which need to be solved for a given background. 
As such, our main results for this part are summarized 
at the end of Sec. IIVI 

There are many gravitational dual theories on the mar- 
ket today. For this reason, it would be useful to have a 
definite prescription which allows for the calculation of 
the bulk spectral function which works for a large class 
of gravity duals. We have taken a step in this direction 
by including multiple scalar fields. Also, as mentioned 
in the introduction, more than one scalar field often ap- 
pears when deriving a five-dimensional effective action 
for a string theory setup. Hence, the methods we have 
developed here could aid future computations in back- 
grounds like those of [271 125] . 

In the second part of the paper, we used the meth- 
ods derived in the first part to compute the bulk spec- 
tral function (at vanishing spatial momentum q) in the 
Chamblin-Reall background. We noticed the intriguing 
fact that in this model the ratio of the bulk to shear spec- 
tral densities is a pure constant. We also derived a sum 
rule for the bulk channel and show that it is different 
from the sum rule of Romatschkc and Son [?D] . The rea- 
son for this difference is not due to asymptotic freedom, 
but rather due to the "hard" nature of the conformal 
symmetry breaking in the Chamblin-Reall background. 
We have also shown that (to leading order in the confor- 
mal symmetry breaking) the left side of the sum rule in 
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Appendix A: Equations of motion 

In this section we write down the relevant background 
Einstein equations and the linearized Einstein equations 
for the perturbations. 



1. Background equations 

There are 3 + n independent background equations, 
corresponding to three metric components gu, fc, g Z z 
and n scalar fields <Iv For later convenience, we write 
these equations as the vanishing of the quantities 



(0) 



with the definitions 



,(o) 



£■ 



(o) 



(0) 







(Al) 



£[ 0) =V{*i,...* n ) 



~-d z (V=gg zz -D L [g xx }lA2) 



AO) 



KK-2?> L [g xx ]V L 



AO) _ 

Co 



AO) _ 




g zz f 
[g X x\ 2 

~gg zz v L [f]) 



d z (V=gg zz K 



(A3) 
(A4) 
(A5) 
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The superscript (0) denotes that these are background 
(zeroth order) equations. The logarithmic derivative T>l 
was introduced in pdj ). A summation over the repeated 
index a is implied. 



2. Linearized equations 

There are 2+n independent linearized Einstein equa- 
tions (first order in the perturbations) corresponding to 
the two metric perturbations A, B and the n scalar per- 
turbations tpi. In |20j . the linearized perturbations for 
this system were given in the more general case of q ^ 0; 
in the general case there are 4 metric perturbations de- 
noted as A, B, C, D. One can get the equations necessary 
for the current work by examining the equations in |20j 
in the limit of q = 0, D = 0, B = C. We will write 
the equations for the Fourier transformed perturbations, 
such as 



A{t, z) = 



dw 
2^ 



A(w, z)e 



(A6) 



with the shorthand A w , B w , tp i w denoting A{w, z), 
B{w, z), ipi(w, z). As in the previous subsection we write 
these equations as the vanishing of the quantities 



with the definitions 



si 1 ' = 



(i) 



(i) 



o 



B,, 



1 



(A7) 



(A8) 



V L [gl x ]A' w + 2,V L [g 2 xx g tt ]B' w - 2$' o¥ / . 



~99 



:dz [V=gg"K] ~ 6w 2 g zz g u B w ,(A9) 



°<S>i — 



1 ~d z [V=gg"<f/ ita ] + ^(A' W + 3B'J 
99 2 



9z 



d 2 V 



w 2 g zz g tt f i ^ w . 



(A10) 



Summation over the repeated index a is implied, and 
we have used the background equations to remove any 
reference to the potential in ( A9 1 



Appendix B: Derivation of on-shell action 

In this section, we provide a calculation of the on-shell 
action in terms of gauge invariant variables, leading to 
the expressions (l22|) and (24). Throughout this section 



we will assume that the metric is written in the coordi- 
nate system (21). Because it is sometimes cumbersome 



to write down square matrices of arbitrary rank, we will 
write any necessary matrices for two scalar fields. It will 



be immediately obvious how to generalize our final result 
for more than two scalar fields. 

To begin, we start with the action ^ and the per- 
turbations defined in (14 1,(15). Next, one expands the 



action to second order in the perturbations. 



1. Bulk term 

We first focus on the "bulk term" : the first term in {[sj) . 
The part quadratic in the perturbations can be written 



5 2 bulk = ^Jd 5 x [Q T MQ + {Q!) T M ZZ Q! + (Q') T M Z Q 
+ (Q"fM 2 zQ + (Q) T M 2t Q + (Q) T M ft Q\ . (Bl) 
Here, Q is the matrix of perturbations 



Q(t,z) 



( A{t,z) \ 
B(t,z) 
<Pi(t,z) 



(B2) 



the dot denotes derivatives with respect to time, and the 
prime denotes derivatives with respect to z. And we have 
introduced several 2 + n by 2 + n matrices:M, M zz , etc. 
We will not write down the components of all these matri- 
ces explicitly here, as only a fraction of these components 
enter into our final result. 

We will work with the Fourier modes, by inserting 

Q(t,z) = J ^Q(w,z)e~* wt . (B3) 

And we will employ the shorthand 

Q w = Q(w,z). (B4) 
The on-shell action becomes 

bulk _ 

4ft 

+ (Q' W ) T M Z Q^ W + QlMjQ'_ w 

+ (Ql) T M 2z Q_ w + QZ M 2zQ-w + QlMcQ-y,} 

with M c defined as 

M c = 2Af + w 2 (M tt + Mj t - M 2t ~ M% t ) . (B5) 

Note that because of the integration over all w, (from 
—oo to oo ) we have taken the opportunity to symmetrize 
the entire expression in w — > — w. Terms which are odd 
in w will drop out under the w integration. Integrating 
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by parts to get an action which depends only on first 
derivatives, one has 



^bulk 



/ dz^ [(Q' w ) T M a Q'_ w + (Q' W ) T MbQ-i 



4k 

QlM'bQ'_ 



QlM c Q. w +d z A(w,z)] (B6) 



with 



2M„ 



M< 



2z 



M a (z) 

M b {z) = M z - 
A(w,z) = (Q'J T M 2z Q 



Ml 



ML 



QlMlQ- 



(B7) 
(B8) 
(B9) 



and we have introduced the result of the integration over 
d 3 x as V3 . One can now employ the Euler-Lagrange equa- 
tions to get the equations of motion for the perturbations. 
In doing so, one should treat Q w and Q- w as separate 
functions. (The reason is that one has in mind a com- 
plex perturbation so that Q- w ~ Q* w [II].) One can 
show that upon application of the equations of motion, 
this entire term in the on-shell action reduces to a to- 
tal derivative term. This procedure has been explained 
many times in the literature (e.g. [21 US] [53] ) hence we 
will not go through all of the details. The result of this 
procedure is 



r>bulk 

<> 2 — 



4k 



dw 
2^ 



(Q' w ) T (M a +M 2z )Q. 
QlM 2 r z Q'_ w + QlM?Q- 



(B10) 



Here z = Zb is a symbolic way of noting that the above 
should be evaluated at the boundary of the space-time. 



2. Gibbons-Hawking term 

We now must add the second term of Q, which is 
itself a pure boundary term. After expanding this term 
to second order in the perturbations, we find 

S§ H = J (Q'fMoHaQ + Q T M GHb Q. (Bll) 

with two new matrices M G n a , M G Hb- Repeating the ar- 
guments of the previous section, by inserting the Fourier 
modes and symmetrizing in w, we find the Gibbons- 
Hawking contribution 



V3 
4k 



dw 
2^ 



(Q 



^ T McHaQ-w 



+ QZ(M GH b + M% Hb )Q- 



QwM GHa Q'_ w 
(B12) 



By an explicit computation, we find that both M 2z and 
McHa are symmetric, and that 



M GHa = -M 2z . 



(B13) 



Indeed, this should come as no surprise, since the pur- 
pose of the Gibbons-Hawking term is to cancel off the 
boundary contribution due to the integration by parts. 



3. Total 

Adding together the bulk and Gibbons-Hawking terms, 
and using (B13l, we find 



S 2 — S 2 

= Yl 
4k 



bulk 



dw 
2^ 



{Q' w ) T M a Q. 



]{Ml 



(B14) 
2M Gm )Q — w ■ 



We are now ready to write down the forms of the matrices 
appearing here explicitly. For the case of two scalar fields, 
the matrices are 4x4, and are given by 



M a 



, 3/2 

J Qxx 






3 








3 


G 














-2 














-2 



and 



Mi 



f 3/2 
1 Qxx 



2M GHb = 
( -3V L [g xx 
dT> L [g xx ] 


V 



(B15) 



(B16) 



3V L {fg : 





3 1 

XX i 



2$; 


-2^\ 




-6$' 2 











J 



It should be evident how to generalize these matrices 
for more than two scalars. Regardless of the number 
of scalars, the 2x2 submatrix in the upper left corner 
always remains the same, since these entries correspond 
to the metric perturbations A and B. For more scalar 
fields, M a simply has more entries of —2 along the 
diagonal. Similarly, (B16) will in general only have 
non-zero entries in the first two rows; one adds entries 



— 2$g... — 2<§>' n to the first row, and 
second row. 



6$' to the 



4. Gauge invariant form 

We would like to write the action involving the gauge- 
invariant variables 2,^. Our strategy is to employ the 
equations of motion to remove all derivatives except for 
Z' m . Specifically, we first use ( A9 1 to remove all instances 



of A', then remove all ip' in favor of Z' m and B' . Finally, 
use (A8| to remove all remaining B' terms. The result is 



So 



V3 f dw g b J x f 
8k J 2?r g' xx 



(i) 



2Z. 



<S>a \ 



-w,z)$' a c (i) 



T^L[g X xf}B- w - V L [g xx ]A- 
-V3 f dw 



f^j (B17) 

gli 2 f {z^ a {~w, z)z' 9a (w, z) + QiiQ-w) 



4k 



2tt 



The terms on the left side vanish on-shell due to the 



equations of motion £| 



0. The term without 
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derivatives (the contact term) contains the matrix £, the 



components of which can be written (in the case of two 
scalar fields): 



-\T> L \f& 

2 
2 



-It>l\/sL] 

k 2 T> L [k 2 gVx ^/J] 



*1 

2 



7] 



2 

k 2 V L [k 2 gt' x 2 

3Ci[o xx 



/o 






' 3V L [g. 
k 2 



AO) 



3V L [g x!c ] £ 2 



(0) 



o o\ 




0/ 



(B18) 



The second matrix here is proportional to the background 
equations of motion, and thus vanishes on-shell. The 
function x{ z ) is defined as 



X(z) 



-t>l [fg: 



3 I 

XX \ 



f 



d z 



fd z 



Qxx 

g'xx 



Aw 2 



+ - 



2„l 

xx 



(B19) 

It is clear that £ is a symmetric matrix (up to terms 
which vanish upon application of the equations of mo- 
tion). For this reason, it does not enter in the compu- 
tation of the spectral density which we consider in this 
work. The matrix £ may be useful in future work, if one 
is interested in computing the real part of the correlation 
functions. Such a computation would also require extra 
counter terms to renormalize the on-shell action. In the 
present work, we do not need to worry about these com- 
plications, as these counter terms are real, and do not 
contribute to the spectral functions. 



The equation (B17l is the final form for the action 



in terms of the gauge invariant variables. Taking the 
imaginary part of the action 

lmS 2 (w) = I(S 2 H-S 2 (-™)), (B20) 

and using the relations for temperature and entropy 

2ng X x(zh) 3/2 



T 



4tt : 



(B21) 
(B22) 



we arrive at the form quoted in the text, Eq. (24) 



Appendix C: Derivation of the equation for for a 
single scalar 

In this section we provide the detailed steps which al- 



low one to derive (46) from (43) in the case of a single 



scalar field. Let us write equation (43) as 



ef = o 



(Cl) 



with 



~-d z [yf^gg zz Z' 9 ] - w 2 g tt Z* 



d 2 V 



-9 



--d, 



gg"&k] . (C2) 



We now add the following combination of the background 
equations ( A2 ) - ( A5 1 to get 



£ (1) + Zj> 
Z V L \g xx 

1 9z 



-d (£ {0) 
k * 



,(0) 



72 9 * 



-9 



the algel 

d z [V^gg zz z^] 



2V L [g xx ]£^' 
(k^g") 2 4 0) ]\=0- (C3) 



Carrying out the algebra, we find that this reduces to 
1 



z„ j^y 4 + jq=9, [V^gg zz fk'} 1 = 0, (C4) 



which is the equation presented in the text ( 46 1 



Appendix D: Chamblin-Reall spectral function - 
comparison with Gubser et al. 

In this section, we compute the spectral density p bulk 
in the Chamblin-Reall model using the method of Gubser 
et al. |14j . The results are consistent with those given 
in Sec. |V A\ though our methods are more generally ap- 
plicable to theories with multiple scalar fields. The no- 
tations and conventions of this section are not the same 
as those given in the rest of the present work. Instead, 
the notation here is the same as that in [2]. At times 
we have inserted some dimensionful factors for clarity. 

The potential is written as 



V($) = Vbe 7 * 



(Dl) 



with < 7 < y/2/3. The background is written in co- 
ordinates where the scalar field is a scaled version of the 
radial coordinate (which we will call r throughout this 
section) : 



ds 2 = e 2A(r) ^_ h{r)dt 2 + rfgB] + £ 



14 



A(r) 

B(r) 

h(r) 
$(r) 



= 1 — exp 



37 2 



6 7 2 

V^Vo(8 - 3 7 2 ) 

G") 



(r - r h ) 



(D3) 
(D4) 

(D5) 
(D6) 



There is a horizon at r = r^. Primes always denote 
derivatives with respect to the radial coordinate. The 
spectral densities are given by 



shear 



hulk 



(h* 12 h' 12 - h^ 2 h 12 ) (D7) 
.^^ yl (K^-hi^). (D8) 



The above quantities are independent of r; the right hand 
side can be evaluated at any value of r. The perturba- 
tions h\ 2 and h\\ are found by solving the linearized Ein- 
stein equations with the incoming wave boundary condi- 
tion. These equations are: 



h'M + h' 12 (r) 



AA'(r) - B'{r) 



h'(r) 
h(r) 



,B(r)-A(r) 

~W) 



Mr) = 0, 



(D9) 



and, 



V ° -4A / (r) + 3B'(r)- h ' ( - 



3A'(r) 



we 



B(r)-A(r) 



h(r) 



Voh'jr) 
6h(r)A'(r) 



h{r 

h'{r)B'{r) 
Kr) 



h n (r) 
(D10) 



These equations hold for any general single scalar gravity 
dual. The key observation is that for the Chamblin-Reall 
background, 



Vo 



3A'(r) 



2B'(r) = 0. 



(Dll) 



With the use of this equality, one sees that the equa- 
tion for h 12 becomes identical to that of hu. Because 
the differential equations are the same, and the bound- 
ary conditions are the same, we conclude that for this 
background only, hu — h\ 2 , and thus, we again see that 
the bulk and shear spectral densities are simply related 



p hulk {w) 
p shear (u>) 



Vo 
A'(r) 2 



9Y 



(D12) 



By comparing the definition of the potentials (48 1, (Dl) 
we see that 



Again we find 



45 

y 



„bulk 



p 



.shear 



128. 



(D13) 



(D14) 



This is in agreement with the result given in the text 
([57). 
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